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ABSTRACT 

The instability of an incompressible three-dimensional boundary layer (that is, one 
with cross-flow) is considered theoretically and computationally in the context of vor- 
tex/wave interactions. Specifically the work centres on two low-amplitude, lower-branch 
Tollmien-Schlichting waves which mutually interact to induce a weak longitudinal vortex 
flow; the vortex motion, in turn, gives rise to significant wave-modulation via wall-shear 
forcing. The characteristic Reynolds number is taken a s a large parameter and, as a con- 
sequence, the waves 5 and the vortex motion are governed primarily by triple-deck theory. 
The nonlinear interaction is captured by a viscous partial- differential system for the vortex 
coupled with a pair of amplitude equations for each wave pressure. Three distinct possibil- 
ities were found to emerge for the nonlinear behaviour of the flow solution downstream - 
an algebraic finite-distance singularity, far-downstream saturation or far- downstream wave- 
decay (leaving pure vortex flow) - depending on the input conditions, the wave angles and 
the size of the cross-flow. 


1. INTRODUCTION 

Our primary concern in this paper on transition is to extend the recently developed 
ideas on nonlinear disturbances in otherwise two-dimensional boundary layers to the more 
general and practical case of three-dimensional boundary layers, i.e. with cross-flow. Al- 
though some linear theory has been written where cross- flow effects have been incorporated 
(such as Stuart 1963, Hall 1986 on predominantly inviscid modes and Stewart and Smith 
1987 on viscous-inviscid modes) little research has been performed in a nonlinear context 
as far as we are aware. Yet the influences of cross-flow are of particular importance in 
numerous aerodynamical configurations, for example on swept wings or near wing-body 
junctions. 

Concerning flows without significant cross-flow, initially at least, in recent years evi- 
dence has accumulated that persistent streamwise vortices can play a key role as an early 
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stage in transition to turbulence of boundary layers and duct flows. In the boundary-layer 
context, the evidence is mainly from experimental observations (Aihara et al. 1965,1969, 
1981,1985, Tani and Sakagami 1962, Klebanoff,Tidstom and Sargent 1962, Bippes and 
Gortler 1972 ),but also from supporting computational work (Wray and Hussaini 1984, 
Spalart and Yang 1987 ). 

The observations have provided the stimulus for recent theoretical studies in this field 
notably for boundary layers (Hall and Smith 1989,1991, Smith and Walton 1989 ) and 
channel flows (Hall and Smith 1988, Bennett, Hall and Smith 1991, Smith and Blennerhas- 
sett 1992 ). The studies are based on the ideas of vortex/ wave interactions: this is the 
situation where two or more wave-like disturbances mutually interact under inertial forces 
in such a way as to induce a longitudinal vortex flow component. This component gen- 
erally has slower streamwise and temporal variations, given that each wave has common 
frequency and streamwise wavenumber. If the vortex also has a significant back-influence 
on the waves, i.e. contributes to their growth or decay, then we have full interaction and 
the waves’ and vortex solutions must be solved together. The interaction may be weak 
or strong depending on whether the oncoming mean flow is altered by a small amount or 
a significant amount respectively. Quite recently, the weakly nonlinear interactions arising 
between a pair of three-dimensional(3D) Tollmien-Schlichting (TS) waves and their in- 
duced streamwise vortices have been studied theoretically for a two-dimensioned boundary 
layer over a flat surface by Hall and Smith (1989). The present work extends the above 
ideas to 3D boundary-layer profiles. Hence our main concern is with the effects of cross- 
flow and the role it plays, if any, in determining the ultimate nonlinear behaviour of the flow 
solution downstream from the input station. This is felt to be a significant problem to 
study mainly because, as mentioned previously, 3D flows are more commonly encountered 
in practice. 

It should be emphasized that nonlinear interactions between vortices and viscous- 
inviscid waves are addressed here with cross-flow present, starting from the interaction with 
zero cross-flow considered in Hall and Smith (1989), Smith and Blennerhassett (1992) and 
then gradually increasing the cross-flow. Interactions involving vortices and predominantly 
inviscid, Rayleigh, waves have still to be extended fully in this way, leading on eventually 
to the study of nonlinear interactions between mean-flow vortices and inviscid cross-flow 
modes. Starts on the latter extension have been made however by Davis (1992) and Brown 
and Smith (1992), based again on the flow structure holding for zero cross-flow as de- 
scribed by Davis (1992), Brown, Brown and Smith (1992), Smith, Brown and Brown (1992). 
One of the main findings, both in the last two works for vortex interactions with inviscid 
waves and in Smith and Blennerhassett (1992) for viscous-inviscid waves, is the existence 
of persistent vortices that emerge downstream of regions of full vortex- wave interaction in 
certain parameter ranges, i.e. pure vortex flow with decayed waves dominates downstream. 
Again, the above is mostly for vortex-wave interactions in the incompressible regime, as in 
the present work. The compressible boundary layer with vortex- wave interaction is treated 
by Bowles, Elliott and Smith (1992), concerning the effects of surface distortions, and again 
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they find persistent vortices to be formed quite commonly as a downstream product of non- 
linear interaction. When that happens, the vortex- wave interaction (upstream) has served 
to alter the mean flow (downstream) to a stable one containing longitudinal vortices. 

An intriguing issue is whether flow properties such as persistent downstream vortices 
continue when cross-flow is added, and, if so, how the parameter ranges are affected. Our 
particular concern is with the effects of cross-flow on vortex/TS interactions. There are 
several kinds of interaction in fact, even within the context of vortices with viscous-inviscid 
(TS) modes. One related kind is discussed by Stewart and Smith (1992), Bowles and 
Smith (1992), regarding flow responses at relatively high frequencies. These are sigificant 
because, among other things, they provide the first theoretical explanation for the Klebanoff 
and Tidstrom (1959) classic path of transition. Indeed, quantitative comparisons between 
the theory and experiments of Klebanoff and Tidstrom (see also Klebanoff, Tidstrom and 
Sargent 1962) in boundary layers and Nishioka et al (1979) in channel flows yield very en- 
couraging agreement, as Stewart and Smith and Smith and Bowles show. The last reference 
also makes comparisons with the strongly nonlinear break-up theory of Smith (1988), for 
a later stage of transition, and again the theoretical-experimental agreement is good in 

quantitative terms. 

As in Hall and Smith(1989),then,we address nonlinear disturbances sufficiently close 
to the first, lower-branch, neutral station. There the TS waves are governed mainly by the 
triple-deck structure, given that the typical Reynolds number, which is defined below, is large 
(Smith 1979). The structure stays intact even with cross-flow present (Stewart and Smith 
1987). Moreover, if the coordinate scales controlling the vortex are taken first to be com- 
parable with those for the wave, the above structure additionally incorporates the induced- 
vortex motion, as we shall see. The nonlinear evolution process is principally contained 
within the lower deck of the triple-deck structure wherein the velocities Uoo(u,t;,iw),the 
Cartesian coordinates ^(x, y,z ), the pressure poo^P and the time loo^t are scaled in 
the form 


[u,v,w,p,x - X 0 ,y,z - z 0 ,t] = [Re 1/8 it,Re 3/8 u,Re 1/s w,Re 1/4 p, (H a - h) 

Re- 3 ^X,Re- 5/6 Y,Re~ 3/6 Z,Re-^ 4 T] 

near a typical 0(1) station x = xq,z — zq .Here looi u oo-> P<xn v oo represent, in turn, the 
typical streamwise length such as the airfoil chord, the flow speed in the outer stream(in 
the x— direction), the fluid density and the kinematic fluid viscosity. The global Reynolds 
number Re = Uoo^oo^ 1 ta ^en to be a large parameter. The oblique TS waves are 
represented by 

E ia = exp[i(aX + /3i, 2 Z/2 — flT)], (1.2<z, b) 

where a,/3i,/3 2 and Q. are real constants. We note that the spanwise wavenumbers /3 i,/ 3 2 
depend essentially on the crossflow evaluated at the edge of the boundary layer. This is 
due to the dispersion relation which stems from the dominant wave motion, as analysed in 
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Section 3. The waves interact nonlinearly to induce streamwise-vortex flow, in the manner 
E^E^ 1 = E 2 , where 


E z = exp[i( fa -f3 2 )Z/2]. (1.2c) 

represents the vortex component. Equally we note the properties E^E^ 1 = E 2 , E 2 E 2 = E\ 
which correspond to the vortex combining with the first wave to provoke the second wave 
and combining with the second wave to provoke the first wave respectively. 

The governing equations are the unsteady, interactive 3D boundary-layer (triple- deck) 
equations, written down in Section 2, and as shown in the analysis that follows that sec- 
tion, these equations cover both the waves’ and the vortex motion, in (1.2a-c), in effect, 
despite the slower scales associated with the vortex. The scales and the flow structure for 
the TS/vortex interaction in a full 3D boundary layer are examined in Sections 3.1 to 3.5 
below, for which a partial- differential system for the vortex flow coupled with an ordinary 
differential equation for each wave pressure is derived. These interaction equations are 
written down in Section 3.6 and special attention is paid to the case of zero starting vortex 
flow,i.e. where the input consists of only the two waves (1.2a, b) superimposed on the 3D 
boundary layer, since then it is possible to deduce an integral form for the downstream 
evolution of the vortex-streamwise wall shear. The interaction equations are addressed 
numerically and analytically in Section 4, and comparisons are noted therein. Three main 
types of nonlinear behaviour are found to occur: an algebraic finite-distance break up; far- 
downstream saturation of the waves’ and the vortex motion; and far-downstream decay of 
the waves’ and the vortex motion. Significantly in the last option it is found that the waves 
die out rapidly in an exponential manner but the vortex decays in a slow algebraic fashion 
and so in a sense we are left with pure vortex motion. The above options depend on the 
input conditions, the wave angles and the magnitude of the cross-flow. It is interesting that 
the numerical results reflect that the first option occurs in the majority of cases. Finally, in 
Section 5, the results obtained in the previous sections are discussed. 

2. THE GOVERNING EQUATIONS 

Substituting the expansions (l.la-h) into the Navier-Stokes equations shows that the 
scaled variables are governed by the unsteady, interactive 3D boundary-layer equations 


ut 4- uux + vuy + w&z = —px + uyy, (2.1a) 

u>t + uwx + vu>y + wwz = —pz + t«yy, (2.1b) 

u x +vy+wz = 0, (2.1c) 

with the boundary conditions 

u = v = w = 0 at V = 0, (2. Id) 

u ~ Y + A, w ~ Y- 1 as Y — » 00 . (2.1e) 
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Here the pressure p and the displacement decrement A are unknown functions of X y Z and 
T y linked via the expression 




(U'd/d( + W'd/9n) 2 Mt,ii,T) 

[( x -() 2 + ( z - v ) 2 }' /2 


d£dr/, 


( 2 . 1 /) 


arising from the main- and upper-deck analyses, where (u e ,0,w e ) denotes the basic flow 
outside the boundary layer. Bars on the integral signs denote the principal value or the 
finite part. We observe that the cross-flow influence on the above system comes exclusively 
from its value at the w upper edge” of the boundary layer and is felt through the pressure- 
displacement law (2. If). No other cross-flow effect is felt within the lower deck because 
the Cartesian coordinates have been chosen in such a way that the direction of maximum 
basic-flow wall shear is in the “x — direction” , implying that the basic flow has no spanwise 
wall shear. 


3. THE ANALYSIS 


3.1 The main scales 

It is well known that lower-branch TS waves (linear or nonlinear) in a 3D boundary 
layer are governed by the triple-deck equations (2.1a-f) cited above (Smith 1979). Here, the 
waves have characteristically small amplitudes in comparison to their fully nonlinear size; in 
particular, if the relative magnitude is h; where h > Re~ m (for any m > 0 ),then a vortex 
flow of 0(h 2 ) is induced through nonlinear wave-coupling. It is found that the vortex- 
spanwise velocity grows logarithmically far from the plate’s surface due to the algebraic 
decay (of 0(Y~ 2 ) ) of the wave-inertial forcing effects there. The singularity is eventually 
damped out in a buffer deck, lying between the main and lower decks, where the shear- 
inertial effects (essentially proportional to , where X v is the modulated streamwise 

length scale over which the vortex and wave-amplitude variations take place) come into 
play. If the relative thickness of the buffer deck is 6(» 1) then convective-diffusive balances 
yield \d% | ~ 6 -3 . There u v ~ h 2 6 3 In ~ h 2 ln£ , where u v ,w v denote the vortex- 
streamwise and -spanwise velocities, and continuity has been used. In Hall and Smith 
(1989) it was assumed (and subsequently confirmed) that sensitive nonlinear interactions 
would happen if the amplitude variations of the wave over the longer scale were controlled 
to some extent by the vortex shear from the buffer deck. Proceeding in a similar vein 
yields the second 6 - X v balance ~ h 2 6 2 ln£, whereupon combining with the first 

such balance gives 

8 ~ h~ 2/5 ,X v ~ h~ 6/s 

where 

r —2 I - 1 / 2 
h = In h j 
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Unlike the related work in Hall and Smith (1989) for zero crossflow, we have not in- 
cluded nonparallelism here, although it is relatively easy to incorporate these additional 
effects as described in Appendix A. 

3.2 The lower deck 

In the lower deck, where y = Re~ 5 / S Y , viscous forces play a prominent role: this is 
readily observed from the unsteady, interactive 3D boundary-layer equations which hold 
here, namely 

Ut + UUx + VUy + WUz = —px + UYYi 
u>T + uwx + viby + wwz = —pz + Wyy, 
ux + Vy + U>z = 0 , 
with p independent of Y, 

and the no-slip condition : u = u = u> = 0aty = 0 

We refrain from applying the outer constraint (2.1e) here because, in effect, the original 
lower deck has split into two separate decks, characterised by the current regime and the 
thicker buffer deck; it is in the latter region where application of (2.1e) takes place. Here 
we express the flow solution in the form 

u = \Y + hLu^ + h*t*\ z (Z)Y + h 2 L 2 u™ + h} l l*LuW + . . . , 
v = hLv™ + h 2 L 2 v™ + fc n/s £s ( i) + . . . , 
w = hLw + h 2 L 2 w^ + h n / s Lw^ + • ■ • , 
p = hLpW + h u / s Lp^ + . . . + h 16 / 5 p( 3 > + . . • , 

f 1 - 1/2 

where L = | — In h > and A Y is the basic-flow shear arising from the main deck. 

The terms superscripted (1) are TS contributions, with the quantities subscripted by ‘a’ 
denoting the second order effects. The induced- vortex contributions have the superscript 
(3) whilst the quantity A3 (Z)Y is the vortex-streamwise shear that arises through feedback 
from the buffer deck. We note that two streamwise length scales are active: the triple-deck 
scale Af, and the modulated scale X = h~ 6 / 5 X, responsible for wave growth or decay 
and vortex variations. Hence dx — *■ dx + h 6 / 5 d x-, effectively. Finally, we expand A as 
1 + h 6 / 5 Aj + • • • , where Ai is real; this is possible provided the lower deck quantities adjust 
to accommodate the local variation of the skin-friction (Smith and Burggraf 1985). 

We now substitute the above expansions into (3.1a-e) to obtain 

i(aY - D)fi ln + v ln = —iapin + uj n yy, (3.3a) 

i(<*Y — D)wi n = — -finpln + W\nYY , (3.3b) 


(3.2a) 
(3.26) 
(3.2c) 
(3. 2d) 


(3.1a) 

(3.1b) 

(3.1c) 

(3.1d) 

(3-le) 
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(3.3c) 


iau in + VinY + -fSnWin = 0; 


i(ay - 0,)u a in + Yu lnX + ia\iYui n + ia j j + w a in 

+ { A*!*!! } + " h) { -kX } = + ““'" yy ’ (3 - 4a) 

i(aY — fi)t5 a ln + Yw ln % + iaXlYwin + | | = — ^nPaln + Wal nYY ,(3.4b) 


IQcUain Uj n j£ ~h V a \ n y “h 2 /^n^aln — 0i 


(3.4c) 


2 

vnuj 2 y + *M*uy + ^33 + -((Sit2iiwJ 2 ~ ^ 2 «i 2 ^n) = w 33 yy , (3.5a) 

ia(ii* 2 u>n - u n t5j 2 ) + tintiij 2 y + s^ny + ^ 1 - /?2 = ™33yy, (3.5b) 

^ 33 y + ~ (^i — ^ 2)^33 = 0; (3.5c) 

for the main TS, forced TS and vortex flow Respectively, where the associated pressures 
Pin>Poin>P 33 are independent of Y, and zero speed is observed at the wall. Here we have 
employed the TS- and vortex- properties described in the Introduction, that is 

u (1) =u u (X,Y)E 1 + u 12 {X,Y)E 2 + c.c., (3.6 a) 

t£> = u a n(X,Y)E i +u al2 (X,Y)E 2 + c.c., (3.66) 

u (3) = u 33 (X,Y)E 3 + c.c., (3.6c) 

etc., and equated coefficients of E\,E 2 f° r the TS equations and the coefficients of E$ for 
the vortex equations. We return to these equations later after we have established some 
necessary results from the buffer, main and upper decks. For now we observe the far-field 
properties 


Uin ~ l,vin ~ Y,u>i n ~ Y _1 ,p ln ~ 1; (3.7a-d) 

Wain ~ 1 j Wain ~ Y,w a ln ~ Y _1 ,p aln ~ 1; (3.8a-d) 

W 33 ~ Y 3 ln Y,t5 3 3 ~ Yin Y,w 3 3 ~ ln Y,p 33 ~ 1; (3.9a-d) 

as Y —+ oo. The logarithmic results in (3.9a-d) rely on neutrality of the waves and on 
imposing a vortex-spanwise-shear constraint at the wall, namely 
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^33y(0) = r F(Y)dY, 

Jo 

(3.10) 

where F is the wave-forcing on the left-hand side of (3.5b). 


3.3 The buffer deck 


This region is needed to adjust the vortex flow to the outer constraint in (2.1e) above. 

We write the velocity and pressure as 


u = h~ 2 / 5 y + /i 4 / 5 (u< 3 > + Aj y) + • • • + hLvS ^ + • • • + h 11 

+ (3.11a) 

v = h 3 ' s LvW + ■■■ + h*/ 5 v^ + . . . + h*l 5 LvW + . . . , 

(3.116) 

<5 = h'l‘La <» + ■ • • + /i 2 u>< 3 > + fc‘V s Iw<» + ..., 

(3.11c) 

p = HLpM + • . • + A 11 ' 5 !#!') + • • ■ + ft“' 5 p< 3 > + ■ ■ ■ , 

(3. lid) 

where (u,v, w,p) satisfy (3.1a-d),but with Y = /i“ 2 / 5 y, where y is 

0(1). 

Thus we generate the sets of equations 


y^x + 

(3.12a) 

y™x } = -Pz\ 

(3.12b) 


(3.12c) 

+ u^) + (* (3) + Ai y)u { x + w (1) (tl ( s 3) + Xi) + 

*2 = 0, (3.13a) 

y(^2 + 1& § ) ) + (* (3) + A iy) t5 x ) = °> 

(3.13b) 

*2 + ) + 4? = °; 

(3.13c) 

yu^ + v (3) = uf) , 

(3.14a) 

- -(3) .(3) 

y w x = w «» 

(3.14b) 

u'x + v 2 + w ( z 3) = 0; 

(3.14c) 


for the main TS,the forced TS and the vortex flows, in turn, where again the corresponding 
pressures are independent of the normal scale. Each set of equations is to be solved subject 
to the external displacement condition (2.1e) where, for each of the above systems, this 
reduces to 
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(fiW.uWtiW) -4 (^(X),^^ J (3 W)> (3.15a-c) 

(u^ 1 ), w^, w^) ~ y~ 3 , (3.15d-f) 

as y — » oo. Here the displacement decrement A has the expansion 

A = h 4 /h A (3) + hLA W + • • • + /I 11 / 5 !! 1 ) + • • • . (3.16) 

For the main TS part we find simple slip-effect solutions hold: 

= A^\ = — yA^x ■, (3.17a,b) 

*»> = { Afa ft + Afeft}+«, P (I >=P <1 >, (3.17c, d) 

l 2 ay 2 ay J 

where 

A (1) = A 11 (X)Ej + A J2 (X)E 2 + c.c. (3.18) 

Consistent lower-deck matching is required as y — ► 0 + , suggesting that 

u (1) — > A^ as Y — ► oo (3.19) 


in the lower deck. 

Next, (3. 13a-c), solved in conjunction with (3. 15b, e), yields in particular 

<*> = i-'-'A"' + A' 1 ', (3.20a) 

vAI = -y(A^> + A^) - (3.20b) 

where 

li 1 ) = A ail (X)E 1 + A al2 (X)E 2 + c.c. (3.21) 

Compatible matching with the lower deck solution will be guaranteed so long as 

u[') _> A 3 A (1) + I ( a 1} as 7 -» oo (3.22) 


in the lower deck. 

The equations for the vortex, (3. 14a-c), are partial-differential being dependent on y 
and X. They are insoluble at present since they depend on the unknown wave-pressure 
terms j>n,pi 2 ? via a slip condition on u>( 3 ) at the buffer-deck wall. Therefore, we must 
determine the dominant wave variations over the longer scale X , and this requires us to 
solve the pressure-displacement law (2. If) for both the leading and forced TS waves, and 
subsequently coupling the results with each corresponding problem in the lower deck. Since 
the vortex pressure is mainly passive here, we need not concern ourselves with the vortex 
motion beyond this layer. 
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3-4 Pressure- displacement relations 


An alternative to solving the Cauchy-Hilbert integral (2. If) directly is to consider 
the Laplacian equation and associated boundary conditions from which the integral was 
derived , namely 

/ d 2 d 2 d 2 \ 

KdX 2 + dy 12 + dZi) P = °’ (3.23a) 

with ^ 

p' — *■ 0 as y' — *• oo, (3.23b) 

p' -* P, Py' {uedx +w e d z ) 2 A as y' — > 0 + . (3.23c, d) 

Here y'(= Re 3 / 6 y) is the scaled upper deck transverse coordinate, and the pressure p' has 
the splitting 

p‘ = hLp'W + • • • + h^^Lp'W + - - - (3.24) 

where 

P ,(1> = Pn(X,y')E, +p ',2 (X,y')E 2 +c.c., 
pj” = p'u(X, y')B, +p’ ll (X,y')B 2 + C.C. 

Insertion of (3.24) into (3.23a-d) gives 

- (a 2 + /3 2 n /4)p' ln = 0, (3.25a) 

with 

Pin — ► 0 as y 1 — » oc, (3.25b) 

Pin * Pin? Piny ' > — (o:U e + f3 n W e / 2) A\ n 3-S ]/ * 0"^*, (3.25c,d) 

for n = 1,2 to leading order. This yields 

(a 2 + /? 2 /4) 1 / 2 p ln (X) = (a + /3 n w e /2) 2 A ln (X), (3.26a, i) 


for n = 1,2 where, without loss of generality, u e = 1. 
The second level equations are 


9V.i 

dy'- 


~ (« 2 + £n/ 4 )Paln = - 2 ia PlnX' 


(3.27a) 
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with 


(3.27b) 

(3.27c) 

as y' —* 0+, (3.27d) 


Pain 0 as y' — > co, 

Pam -* Pain as y' -> 0 + , 

Palny 1 ~ ( a + PnW e /2) 2 A in + 2i(a + f3 n W e /2)A in x , 

for n = 1,2 which we note is absent of any vortex forcing and therefore linear. In fact only 
the slow stream wise modulation of the main waves serves to drive the secondary waves. 
The solution of (3.27a-d) is 

(a 2 + ft /4) 1 / 2 p aln (X) = (a + P n w e /2) 2 A a i n (X) - 2 i(a + p n w e /2)A lnjt 

ia(a + p n w e /2) 2 - ( 3 . 28 a, b) 

+ WTpi! 4 ) ln *’ 


for n = 1,2. 


5.5 Neutral eigenrelations and pressure-amplitude equations 

Having established the necessary upper-deck displacement laws, we now return to the 
lower deck and solve the leading TS and forced TS systems of equations therein, together 
with (3. 26a, b), (3. 28a, b) above. So for the leading waves we have 

i(aY — fl)ui n + v ln = —iapin + tti n yy, (3.29a) 

i(aY - fl)u)i n = - ^PlnPln + WlnYY, (3.29b) 

£ 

lO'Ujri + VlnY 4" ^ PnWln = Oj (3.29c) 

with 

Uln = «ln = Win = 0 at Y — 0 (3.29d) 

Win ~ iiln, Win ~ Y~ 1 as Y -» oo (3.29e) 

and 

(a 2 + = (a + /3 n w,/2)*A ln , (3.29/) 

for n = 1,2. From this we can deduce the dispersion relations 

= (i a )V*(a* + /3ll4)'l\a + 0„W'/2)\ (3.30a, b) 

K 

for n = 1,2, where 
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i 1 / 3 ^ 


£o = 


a 


2/3 ’ 



and Ai is the Airy function. The neutrality of the waves implies the values 


k d 2 

where di « 2.3, d 2 % 2.3, (Lin 1955, Smith 1979,Drazin and Reid 1981), so that (3.30a, b) 
become 

« 5/3 (<*iM) = (a* +0l/4) ,/2 (c. + l3 n W'/2)\ (3.310,6) 

for n = 1,2. Notably the presence of cross-flow precludes any symmetry in the a - (3 
plane e. /5j = — j3 2 is impossible for any streamwise wavenumber unlike the instance of 
zero cross-flow where perfect symmetry occurs in the a— axis. 

For the forced TS waves the governing equations are 


i(aY i&XjYuin + ia ^ ^ + Vain + Vln^l 

+ { ^ 33^11 } + 2^ } = ~ ta P*in -PmX + “ai nry, (3.32a) 

i(aY - Q)w aln + Yw lnX + iaXiYwm + j J = -~/3 n p aln + u> alnyy , (3.32b) 

iau aln + u lnX + v alnY + ^(3 n Wa in = 0, (3.32c) 


with 


“aln ~ 


^aln — ^aln — V^aln — 0 


at Y = 0, 


f -^ 33^12 
\ ^ 33-^11 


} 


"b A a in,W aln ~ Y 


as Y 


oo, 


(3.32d) 

(3.32e) 


and 

(« 2 +0»/4) ,/2 p.i» = (a+(3 n w e /2 ) 2 A aln (X) - 2i(a+p n w e /2)A lnX 

ia(a + (3 n w e /2) 2 - (3.32/) 

(« 2 +/3 2 /4) 

for n = 1,2. After much working, we deduce the pressure-amplitude equations 
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a \ jjr + ^i^lPn + C1X33P12 = 0 , 

a 2 + C 2 A 33 P 11 = 0 , 


(3.33a) 

(3.33b) 


hold, where 


= 2 Bnrj|or L P - b _ 1/2| 


-2 


3a 3 A 2 / 3 


; _ 2B n r 2 j () r- l D 
n 3aA 5 / 3 


r n J « 2 

1 

tol 01 
+ 

tol ^ 

1 - . — 1 

(3.34a) 

5Bl /2 

Za ’ 


(3.34b) 


2(^ 1 Jg 2 ) 1/2 T^r 1 ^ (TnTn(3lf3 2 r 2 \ 

3 a 2 r 2 r 2 8 a 2 / 

-1 1/2 /5t^(5iS 2 ) 1/2 Z(Tnf3i02 \ 

■“ S " ( Irfrl + 8 > 

for n = 1,2. Here B n = (a 2 + fil/4),T n = (a + f3 n w e /2), a n = (1 - /3i/3 2 //3 3 ), D = 
1 + *6)A4i(£o)j A = za,ri = Az(£o)Mf f (£o )> r 2 = and £ 0 ?^ are as defined above. 

The pressure-amplitude equations illustrate that the growth rate of each wave is af- 
fected by the basic-flow shear correction (Ai) and the buffer-deck vortex- streamwise shear 
(A 33 ); the cross-flow influence is reflected by the interaction coefficients We 

are prevented from solving the amplitude equations as they stand because the vortex (and 
hence A 33 ) has implicit dependence on the wave pressures pu ,pi 2 , as mentioned in Section 
3.3 above. Instead we must solve the wave- and vortex-equations interactively, although 
firstly we need to ascertain the boundary conditions for the vortex at the buffer-deck wall. 

We know from above that, in the lower deck, wave-inertial forcing provokes logarithmic 
growth in the vortex-spanwise velocity component. A more detailed evaluation comes from 
substituting the asymptotic properties 



c' n =W(o&-V>( 


Vln = 


= —(iaAi n )Y + iflAin — ia(l + d” 




(3.35a) 

(3.35b) 

(3.35c) 


of the leading waves into the lower-deck spanwise- vortex equation (3.5b) above. It is found 
that 


where 


™33 


~ — iKpiyj>\ 2 InY, for Y 1, 

* = 5 <A-A)(i + ££). 


(3.36) 

(3.37) 
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The logarithmic growth in W 33 is known to induce the other properties 633 ~ Y 3 In Y, v 33 ~ 
Yin Y as Y — ► 00 , in the lower deck. Therefore, the inner constraints for the vortex in the 
buffer deck are 

U 33 — 1 ► 0,633 — 1 ► 0 , 1633 — > — iKpiip1 2 as y — ► 0 + . (3.38a — c) 


3.6 The interaction equations 


In summary, the nonlinear vortex/TS interaction is embodied in the equations 



yU x + V = u„. 

(3.39a) 


yW X = Wyy, 

(3.39b) 


Ux + Vy+ l -( 0 1 -( 3 2 )W = O, 

(3.39c) 

with 

U(X,oo) = A, W(X,oo) = 0 , 

(3.39d) 


C/(X,0) = 0 , W(X, 0) = -iKP 11 P* 2 , 

(3.39e) 

and 

+ blXlPll + C 1 A 33 P 12 = 0 , 

(3.39f) 


dPi 2 

+ b 2 XjPl2 + C2-^33^11 = 0 , 

(3.39g) 

where -0,33,^33, 
bn — b' n /a' n ,c n 
simplify to 

W33, A 33 , Pn,Pi2 have, in turn, been replaced by U, V, W, A, Pu, P 12 , and 
= c' n /a' n (n = 1 , 2 ). Defining r = Uy (so that r(J?,0) = A 33 (X) ),( 3 . 39 a-c) 


T yy ~ V T X = P2 )W, 

(3.40a) 


Wyy ~ WX = 0 , 

(3.40b) 

with 

t(X, 00) = W{X, 00) = 0 , 

(3.40c) 


rs(X,0) = 0, IY(J?, 0 ) = -iKP 11 P 1 \. 

(3.40d) 


In principle, we may solve (3.40a-d) collectively with (3.39f,g),given some prescribed input 
conditions in X. Analogous equations to (3.39a-g),(3.40a-d) were obtained by Smith and 
Blennerhassett (1992) for zero cross-flow, where the authors corrected the original zero 
cross-flow “interaction equations” in Hall and Smith (1989) . In both papers, a partial- 
differential finite-scheme was applied directly to the vortex-wave equations. 
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It is possible to reduce the r — W system above to an integral equation for A33, 
however, if we consider the special case of zero-input vortex flow, i.e. W = r = 0 at X = 0 
(without loss fo generality). We apply the Laplace transform in X to (3.40a-d) and obtain 


W = - 


%KP 

^4i(0) 


Ai(s 1 / 3 y), 


(3.41a) 


and 


where 


- KP (ft 


2Ai(0) 


(3 2 )s 2/3 Ai(s 1/3 y), 


(3.416) 


1*00 /»oo 

'■(s,y) = / T(X,y)e- a *dX , W(s,y) = W(X ,y)e~' X dX , 
t/o t/o 


P(s 


>-Jf 


P 11 (X)P 1 * 2 (X)e- i;f dX 


are the Laplacian transforms of r, W, Pi 1 P12 respectively. Hence, (3.41b) evaluated at 
y = 0 and inverted gives 


r* 

A 33 (X) = MKfi / P xl {u)P; 2 (u){X -u)-'! 3 du, 
Jo 


(3.42) 


where ft = (/?i — /?2)/2,M ~ 0.54. Then unifying (3.42) and (3.39f,g),and prescribing 
values for Pu,Pi2 at = 0, we can determine the flow solution for X > 0. Numerically, 
our task is much easier, since we have eliminated one variable ( y ) entirely and we do 
not therefore need to resort to the potentially difficult and computationally expensive 
two-variable finite-difference schemes. On the other hand, for non-zero vortex input the 
same Laplace transform scheme for (3.40a-d) would yield additional terms, generally triple 
integrals, on the right-hand-side of (3. 42), and in this case the finite-difference scheme would 
possibly be the better choice. 

For all the ensuing weakly nonlinear analysis, we concentrate on the case of zero- 
input vortices. Computations have been performed for a full 3D boundary layer, i.e. 
w € = 0(1), and these are presented in the next section. 


4, COMPUTATIONAL SOLUTIONS 


We applied a predictor-corrector scheme of second-order accuracy to advance the wave 
pressures in distance and a trapezoidal rule to calculate the vortex-shear integral at each 
station. This procedure proved to be stable and accurate for suitably small step lengths. 
Interaction results were obtained for sample values of a and w e and all starting at X = 0 
upstream of the neutral TS point. The input value for each wave pressure was fixed at 
0.1. The basic-flow correction Ai was taken as positive and therefore normalised to +1. 
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Three main types of nonlinear behaviour are found to occur downstream and, in light 
of the numerical findings, we now address these analytically. Firstly, there is the situation 
(figures 2,5,8 and 9) where the flow solution develops an algebraic singularity at a finite 
position downsream,say as X — » X, ; the orders of magnitude suggest the scalings 


j Pn~r 5/ V* l(c) |P„| + ---, (4.1a) 

P 12 ~C 5/ V< mc) |Pi 2 | + ---, (4.1b) 

■^33 ~ H , (4.1c) 

where £ = ( X 3 — X ), and 0 < ( < 1, Here the real-valued phase factors <f>i,<j >2 expand as 

~ <^io In £ + <f>\\ + <^i2C + • • • > (4.2a) 

<f>2 ~ 4 > 20 In C + 4> 21 + ^ 22 C + • • ■ 5 (4.2b) 


where <f>i 0 , <^ n , <f> 12 , (f> 2 o, <£ 21 , ^> 22 , • • • are unknown constants. The analysis produces the 
solvability conditions 


K(c lr - 6 0 c u ) < 0, (4.3a) 

K{c 2r - S 0 c 2 i) < 0, (4.3b) 

where c lr , cn denote the real and imaginary parts of c\ in turn and likewise for the other 
interaction coefficients. Also £o(== <^io ~ <^20 )> the dominant phase difference of the pres- 
sures, satisfies the cubic equation 


a 3^o a 2^o a i$o + a o — 0, 


(4.4) 


where 


a 3 = 6ciiC 2 i, a 2 = ll(ciic 2r — ci r c 2 j), 
a i = — 2(8ci r c 2 r + ciic 2 j), ao = 5(ci r c 2 i — ciic 2r ). 


(4.5a,b) 

(4.5c,d) 


form 


Secondly, there is the option of saturation where the nonlinear quantities develop the 


P11 ~ |Pn|e ieiJf , P12 ~ |Pi 2 | e * , 

as X — > 00 (see figures 6 and 7). The analysis yields the criteria 

-^[e-^d] < 0, ^-3?[e iiry,3 c 2 ] < 0, 

b \r 02r 


JO 7 X 


and 


A J- c. a[e-”/ »c] 


- (&1« — i)) > 0, 


(4.6a, b) 


(4.7a, 6) 


(4.7c) 
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where £<>(= $1 — # 2 ) 1S the wave-pressure phase difference at the point of saturation, and 
S denote the real and imaginary parts of their enclosed quantities, respectively. 

Finally, there is the possibility that both waves will decay far downstream in an expo- 
nential manner under linear forces, i.e. acquire the form 

P 11 ~ |Ai|e~ kljt , P 12 ~ \Pi 2 \e~ b:,X (4.8a, b) 

as X — * 00 . This is clearly valid so long as 

b lr > 0, b 2r > 0. (4.9a, 6) 

The vortex-shear decays in a much slower (algebraic) fashion because of its dependence on 
the history of the flow solution. Figures 3 and 4 illustrate this option. 

The values of the wave angles and the cross-flow upon which the interaction coefficients 
depend generally determine which option will come into play downstream but no explicit 
constraints on these parameters have been deduced due to the complexity of the analysis. 
The numerical results do indicate however that an algebraic singularity only occurs if one 
of the modes is the “lower” branch cross- flow mode. (See Appendix B for a description 
of the modes and a summary of the linear neutral stability in the presence of cross-flow.) 

5. FINAL COMMENTS 

It has been shown in the preceding sections that significant interactions between two 
TS waves and an induced or incident vortex in a cross-flow boundary layer can lead to three 
types of behaviour downstream. An interesting point is that the results are very much in 
line with the related work on channel flows and boundary layers for zero cross-flow by Smith 
and Blennerhassett (1992). Firstly, the flow solution breaks down in an algebraic fashion at 
some finite position downstream. This points to the entry of a stronger and probably fully 
nonlinear stage more locally. Secondly, the flow solution may saturate far downstream, i.e. 
the flow quantities may asymptote to finite values as the streamwise coordinate becomes 
large. The computations show however that this second option is fairly infrequent and only 
seems to occur in a relatively small range of angles for each given cross-flow. This is consid- 
erably different from the Smith and Blennerhassett results, where most of the wave angles 
considered lead to the saturation case. Moreover, in their case the wave pressures decay 
to zero during saturation whereas the wave pressures in our case asymptote to non-zero 
values. Thirdly, the TS waves may decay exponentially at downstream infinity, essentially 
under linear forces, leaving pure vortex flow. (We note,however,that the magnitude of the 
vortex flow is smaller than in the previous region of full vortex/wave interaction, because in 
that region it decayed algebraically.) The numerical results tend to indicate that this third 
option happens only when both of the TS wave angles are greater than — arctan^jT 1 ) 
to the x — direction (which we recall is the direction of the surface shear stress) . This, in 
turn, is possible only when the four-mode criterion, i.e. w e < 0.20 approx., as discussed in 


17 



Appendix B,is in effect, and neither mode is the “lower” branch cross-flow mode, whose an- 
gle is always less than - arctan(tt?7 1 ) to the x- direction. In contrast, the first and second 
options seem to occur only when one of the modes is the “lower” branch cross-flow mode. 

The neutral stability curve that stems from linear TS theory has been illustrated for 
the cases w e = 0, 0 < w e <C 1 and w e = 0(1) in Appendix B. It can be seen that for 

all non zero cross-flow values the curve is open and unbounded, reflecting the property that 
the chosen wavenumber scales are not sufficiently large to capture the complete neutral 
curve, i.e. the “closure” of the curve. There the relevant streamwise scale in particular is 
0(J?e“ 3 / 7 ) instead of the triple-deck scale 0(Re ~ 3 / 8 ), and corresponds to having signifi- 
cant vertical wave- acceleration entering the boundary layer and inducing a wave-pressure 
“jump” across the extremes of the boundary layer. If the jump is too large then the bound- 
ary layer cannot support neutral solutions (see Davis 1992). It is hoped that the weakly 
nonlinear theory will be extended to that regime in the future. 

The current regime of triple-deck wavenumbers is of further interest because it enables 
comparisons to be made between the linear natures of the two cross-flow modes and the 
two regular 2D-type modes. More importantly, for cross-flow values less than 0.20, these 
modes co-exist, enabling us to consider the influence of each of the six possible pairs of 
modes on the weakly nonlinear interaction. In addition to the 0(1) computations pre- 
sented above, some computations were also performed for the instance where the cross-flow 
parameter w e is asymptotically small, based on the corresponding analysis in Davis (1992). 
In that analysis, it was shown that the two cross-flow modes expand in the form 


(3 ~ (~2a)w e 1 ±$iw e 1/2 + 0(1) asw e ->0 + , (5.1a, b) 

where 

Pi = 2(d 1 /d 2 ) 1/2 a 1 ^, (5.2) 

illustrating that these modes are nearly normal to the x— direction. Also it can be seen 
that the cross-flow influence remains despite the smallness of the cross-flow magnitude; 
this contrasts with the two regular modes which have the approximate 2D- boundary-layer 
forms 

P ~ ±/?0 + Pi w e + 0{w\) as w e -* 0+. (5.3a, b) 


Here /? 0 (> 0) is the “zero cross-flow” mode of Hall and Smith (1989) which satisfies the 
neutral equation 


d j 

d 2 


= a 1 ' 3 (a 2 




(5.4) 


and 




(5.5) 


The computations for small cross-flow showed that the option of exponential decay would 
occur in the weakly nonlinear interactions unless one of the modes was the “lower” branch 
cross-flow mode. In the latter case the flow solution would blow up algebraically. These 
results are also largely consistent with the computations for 0(1) cross-flows. 
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We conclude therefore that for cross-flow values less than 0.20 the nonlinear TS/ vortex 
interaction either undergoes an algebraic finite-distance blow-up or else the exponential- 
wave-decay /algebreiic- vortex-decay option comes into play, depending entirely on whether 
or not the “lower” branch cross-flow mode is one of the two modes involved in the interac- 
tion^ discussed previously. For cross-flow values exceeding 0.20 however, the flow solution 
breaks up algebraically in all cases, except possibly in a small angle range where saturation 
occurs. 

On a final matter, it is clear that this and probably the other vortex/TS-wave in- 
teraction structures possible (see also Section 1) tend to stay essentially unaltered by the 
additional presence of 0(1) cross-flows, i.e. where the edge velocities u € ,w e are comparable. 
By contrast, recent work by Davis (1992), Brown and Smith (1992) on the vortex/inviscid- 
wave interactions shows that even small amounts of cross-flow can substantially alter the 
original flow structures set up by Hall and Smith (1991) (for long scales), Smith, Brown and 
Brown (1992) (for short scales) for such interactions with zero cross-flow. There appears 
to be much to be clarified there. 
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APPENDIX A. EFFECTS OF NONPARALLELISM 
ON THE TS/ VORTEX STRUCTURE 


The global basic-flow velocity ( u,Re 1 ^ 2 v,w) has the expansion 


u = u(y) + (x - x 0 )u b (y ) H , 

Mi) 

v = Vb(y) H , 

M2) 

w = w(y) + (x - xo)wb(y) H , 

M3) 


as x — ► in the main deck. Nonparallelism becomes significant when the second terms 

in u,w and the first in v, start to drive the vortex-induced waves over the longer scale X. 
Hence, we seek the balance 


Re~ 3/6 h- 6/5 ~ h 6/5 , i.e. h ~ Re~ 5 / 32 . 

The subsequent alteration to the interaction equations is the insertion of additional terms 
(&iAbX)Pii and ( bi\bX)P \2 in (3.39f,g) respectively, where A b = U(, s (0). This may possi- 
bly lead to a substantial change in the development of the flow solution, especially if the 
nonlinear structure avoids breaking up after a finite distance downstream, but on the other 
hand nonparallelism may be passive. Much depends on the wave angles, the size of the 
cross-flow and the input conditions. 

We observe that for h Re~ 5 / 32 , non parallelism is negligible. Thus with hind- 
sight, the analysis of Section 3 is seen to be valid in the regime 

Re~ 5/32 < h < 1. 
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APPENDIX B. EFFECTS OF CROSS-FLOW 
ON THE LINEAR TS STABILITY 

The neutral stability curve that corresponds to the linear 3D dispersion relation 

is illustrated in Figure 1 for the cases w e = 0(1), 0 < w e <C l,w e — 0. 

In particular, there are some striking differences between the last two curves: firstly 
the zero cross-flow curve is symmetric and extends across the entire f3— range, whereas 
the small cross-flow curve remains asymmetric and has a finite cut-off value for some 
/3 > 0, above which no neutral solutions exist; secondly there is a vast difference between 
the respective maximum a— values. When w e = 0 this value is (dj/^) 3 / 4 ~ 1 and occurs 
when (3 = 0 (i.e. for 2D waves) but for all non- zero values of the cross-flow it is removed, in 
effect, to /? = — oo and becomes infinite (pointing to the existence of important long-scale 
instabilities). Obviously the curves do not match uniformly (even though the regular 
“upper” branches coalesce) and an intermediate matching regime exists. Another feature 
of the small cross-flow curve distinct from the 2D curve is the presence of two “extra” 
modes, these being purely cross-flow generated. This points to the possibility of having 
six pairs of interactions for any given ex lying inside the critical interval (0,1), (where we 
have given {d 1 /d 2 ) the appromixate value of 1). This interval becomes thinner as cross- 
flow increases until eventually it vanishes for a critical cross-flow value, above which only 
the cross-flow modes exist. (This critical value has been determined to be (-v/24) -1 (~ 
0.20), where the “upper” branch has no maxima or minima.) Asymptotic solutions have 
been obtained for the weakly nonlinear interactions in the limit of small cross-flow (Davis 
1992) . 
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(a) TS pressure moduli (|Pn|, P 12 1) versus X. 



(b) Vortex shear modulus IX33I versus X 

Figure 2. — Nonlinear-interaction computed results. Wave angles at 
86.42°, -89.59°. Cross-flow at 0.01, grid AX = 1 0- 5 . Initial TS 
pressure moduli both 0.1 , zero-input vortex flow. 



0 5 10 15 20 

X 

(b) IX33I versus X. 


Figure 3. — Nonlinear-interaction computed results. Wave angles 
8-j = 65.56°, 8 2 = -67.80°. Cross-flow is 0.01 , grid AX = 5 X 10“ 3 . 
Initial |P1l|,|Pl2| both 0.1, |x 33 | = 0. 






(a)|Pn|.Pi2l versus X. 

2 



0 

<a)|Pn|.|Pi2| versus X. 


(t>) |X 33 | versus X. 

Figure 4. — Computational solutions of the nonlinear interaction. 
Wave angles 6 1 = 37.72°, 0 2 = “65.36°. Cross-flow is 0.1 , grid 
AX = 5 X 1 0-3. Start |P1 1 1 , |P1 2| both 0.1 , |\ 33 |= 0. 


(b) IX33I versus X. 

Figure 5.— Computational results of the nonlinear interaction. 
Wave angles at 6 1 = -76.31 °, 6 2 , -85.70°. Cross-flow is 0.1 , 
grid AX = 1 0“ 4 . Initial |P1l|,|P12|both 0.1, Ix^l = 0. 
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